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Abstract 

We present an "algebraic treatment" of the analytical Bethe Ansatz. For this purpose, we 
introduce abstract monodromy and transfer matrices which provide an algebraic framework for 
the analytical Bethe Ansatz. It allows us to deal with a generic gl(M)-spm chain possessing on 
each site an arbitrary g£(AA)-representation. For open spin chains, we use the classification of 
the reflection matrices to treat all the diagonal boundary cases. 

As a result, we obtain the Bethe equations in their full generality for closed and open spin chains. 
The classifications of finite dimensional irreducible representations for the Yangian (closed spin 
chains) and for the reflection algebras (open spin chains) are directly linked to the calculation 
of the transfer matrix eigenvalues. 

As examples, we recover the usual closed and open spin chains, we treat the alternating spin 
chains and the closed spin chain with impurity. 
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1 Introduction 



The investigation of the integrable quantum spin chains was initiated by H. Bethe in 1931 [1] where 
he studied the closed spin 1/2 Heisenberg chain [2]. Since then, numerous generalisations of this spin 
chain have been introduced: anisotropic XXZ spin chain [3-5]; spin 1 chains [6,7]; alternating spin 
chains [8,9]; spin chains with higher spins [5,10-15]; spin 1/2 chains with spin 1 impurities [16-18]. 
Correlation functions of this type of spin chains have been computed in e.g. [19-25]. The framework 
of integrable open spin chains has been developed in [26-29]. 

There exist different motivations to study generalisations of the integrable spin chains. First, 
they describe dynamics which can be computed exactly of quantum mechanical models. Indeed, new 
models have been investigated to describe theoretically crystalline material in order to compare with 
the experimental data (e.g. for the crystals MnCu ( S 2 C 2 2 ) 2 {H 2 0) 3 see [30], (VO) 2 P 2 7 [31-35] or 
Cu 2 (OH) 2 C03 [36]). Spin chains allows one to treat some limits of other models as the Hubbard 
model [37], the quantum chromodynamics theory [38] or integrable relativistic quantum field theories 
[39-42]. Finally, recent developments in AdS/CFT correspondence have also put spin chain models 
in foreground of string theory [43-45]. 

The increasing number of applications urge us to seek for a complete treatment of (closed and 
open) spin chain models. Different schemes exist for dealing with these problems, most of them re- 
lying on the Bethe Ansatz (coordinate, algebraic, analytical or thermodynamical for the main ones). 
We present here a formulation of the analytical Bethe Ansatz for gl(N) (closed and open) spin chains 
whatever the representation at each site the quantum spins belong to. In particular, we unify by this 
way all the generalisations of the XXX model. 

The main results of this paper are the following: 

- The determination of the Bethe equations for a closed spin chain model where each quantum 
spin is represented in an arbitrary representation of gl(Af) (called generic closed spin chain). 

- The computation of Bethe equations for any open spin chain constructed from an arbitrary 
diagonal reflection matrix and a generic closed spin chain. 

- For each of the above mentioned models, the calculation of the underlying symmetry and the 
integrability of the models. 

This paper consists of two main sections. The first one is devoted to the study of closed spin chain, 
while the second one deals with the case of boundaries. The structure of these two parts is similar. We 
recall first the algebraic settings (Yangians or boundary algebras) for the monodromy matrix. Then, 
we use the classification of the representations of these algebras to compute a represented transfer 
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matrix. Next, we use generalisation of analytical Bethe Ansatz to obtain the Bethe equations. 
Finally, we work out some examples. 



2 Closed spin chain 
2.1 The R matrix 

We will consider the gl{H) invariant R matrices [46,47] 



R ab {\) = lN®lN- f ^y L , (2.1) 
where V ab is the permutation operator 

N 

V ab = J2 E ij ® E a ( 2 - 2 ) 

and h is the deformation parameter. It is usually set to 1 in the context of quantum groups (Yangians), 
and to —i when dealing with spin chain models: here, we leave it free. are the elementary matrices 
with 1 in position and elsewhere. From the algebraic point of view, the value of h is irrelevant 
(provided non- vanishing) . It is in general set to 1 when studying Yangians, while it is set to — i in 
the spin chains context. Here, we leave it free. 
This R matrix satisfies the following properties 

(i) Yang-Baxter equation [4,5,46-48] 

Rabija — Aft) i? ac (A a ) R bc (X b ) = Rbc(^b) Rac(^a) Rab(^a — Afe) (2.3) 

(ii) Unitarity 

Rab(X) Rba(-X) = C(A) In ® Ijv , (2.4) 
where R ba (\) = V ah R ab {\)V ab = i?^ b (A) = R ab (X) and 

™-H)H)- m 

It obeys [A a A b , R ab (X)} = for A E End(C M ). 

The R matrix can be interpreted physically as a scattering matrix [3,4,49] describing the interac- 
tion between two solitons (viewed in this framework as low level excited states in a thermodynamical 
limit of a spin chain) that carry the fundamental representation of gl(Af). 
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2.2 Yangian y(gl(J\f)) 

We present in this section some definitions and properties of the Yangian [50] associated to the Lie 
algebra gl{M) that will be used in the following. 

The Yangian y(gl(J\f)) is the complex associative unital algebra with the generators {T^ |1 < 
i,j < A/", n G Z> } subject to the defining relations 

. T m] ~ PS , Tt 1] ] = T« TW - r W r« , (2.6) 

where r,sG Z> and T^- = (5^-. 

The i? matrix previously introduced allows us to encode the Yangian defining relations in a simple 
equation, called FRT exchange relation [51] 

R a b(K — A&) T a {\ a ) T b (\ b ) = T b (X b ) T a (X a ) Rabija — A?,) , (2.7) 

where the generators are gathered in the following matrix (belonging to End(C^) ^)y(gl(Af))[[\^ 1 ]]) 

T(A) = E * ® T ^ A ) = E ^ ® E J ^ = E ^ rW • (2-8) 

i,j=l i , jr = 1 r>0 r>0 

Using the commutation relations (|2.7|) . it is easy to show that generates a gl{M) algebra. 

In order to construct representations of y(gl(Af)), the following algebra homomorphism from y(gl(J\f)) 

to U(gl(J\f)) (universal enveloping algebra of gl(N)) will be used 2 

T^X) ^ 5 {j - ^ , (2.9) 

where {e^} is a basis of the Lie algebra gl(N). The Yangian of gl(N) is a Hopf algebra with the 
coproduct given by 

A:y( g i(Af)) — :%W)®W(-A0) 

fe=l 

The coproduct is the cornerstone to deal with the tensor product of representations. We define also 
by recursion A^ = (A ® id ® n_2 )A (n_1) for n > 2 and A( 2 ) = A. 

The quantum determinant qdetT(X) is a formal series in A -1 with coefficients in y(gl(N)) defined 
as follows 

qdetT(X) = s 9n(<?) Zv (1) (A - W + h) ■ ■■T NMM) {\) , (2.11) 

2 To be compatible with the pseudo-vacuum as usually defined in the study of spin chain models, the convention 
used here for the homomorphism differs from the one introduced in [52]. The link between the two conventions is 
provided by the Yangian automorphism T(A) i — ► T (— A), where * is the usual transposition. 
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where &n is the permutation group of Af indices. A well-known result (see e.g. [53]) establishes 
that the coefficients of qdet T(A) are algebraically independent and generate the centre of y(gl(Af)). 
It is important for the following to realise that the quantum determinant represented in any finite- 
dimensional irreducible representation will be proportional to the identity matrix. 

There exists an equivalent definition of the quantum determinant which will be used in the 
following as well. Let A m be the antisymmetriser operator in (C )® m , i.e. 

A m (e n ® ■ ■ ■ <g> e lm ) = i s 9 n ( a ) e v (1) ® • • • ® e ia(m) , (2.12) 

where {ej|l < i < Af} is the canonical basis of and 1 < i\, . . . , i m < Af. The antisymmetriser is 
a projector in (C^)®" 1 . It has the remarkable property: 

Proposition 2.1 [54] The following identities hold 

A m T l {X)---T m {\-mh+fi)A rn = T m (\-mh+h)---T 1 (\)A m (2.13) 

= A m T l (X)---T m (\-mh + h) . (2.14) 

When m = Af, the antisymmetriser becomes a one- dimensional projector and one can show [54] 

qdet T(A) A N = T M (\ -hAf + ft) ■■■ 71(A) A M . (2.15) 

The relation ()2.15|) can be used as an equivalent definition of the quantum determinant. 
To study spin chains, we will use the following automorphisms of y(gl(Af)) 

(i) Inversion 

inv : T(A) i — ► T-\-\) (2.16) 

(ii) Shift 

s a : T(A) i — > T(A + a) , a G C . (2.17) 

One can compute the elements of T _1 (A) in terms of T(A) using the following formula 

T-\X -hAf+h) = {qdet T(A)) -1 T*(A) , (2.18) 

where T*(A) is the quantum comatrix, i.e. its entries Tj*(A) are (— l) l+J times the quantum determi- 
nants of the submatrices of T(A) obtained by removing the i th column and j th row. 

2.3 Algebraic transfer matrix 

In the following, in order to construct spin chains, it will be necessary to deal with the tensor product 
of £ copies of the Yangian. For 1 < i < £, we denote by £ ai (A) G End{C M ) <g> y{gl{Af)) one copy 
of the Yangian which acts non trivially on the i th space only. The space a, always isomorphic to 
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End(C^) in the present paper, is called auxiliary space whereas the space i is called quantum space. 
Obviously, C ai {\) satisfies the defining relations of the Yangian 

Rabija — Aft) C ai (X a ) C bi (X b ) = £bi(X b ) £ai(X a ) Rabija ~ X b ) . (2-19) 

Let us stress that the matrix C a i(X) is local, i.e. it contains only the i th copy of the Yangian. On 
the contrary, thanks to the coproduct, one constructs a non-local algebraic object, the monodromy 
matrix 

r.(A) = A^(£(A)) = £ al (A) C a2 (X) . . . C a£ (X) e End{C«) ® {y{gl{N))) m . (2.20) 

Let us remark that the quantum spaces are omitted in the LHS of (|2.2()|) . as usual in the notation of 
the monodromy matrix. The entries of the monodromy matrix T a (X) are given by 

M 

T *M) = Yl L ^ ® L kik 2 (X) ® ... ® L k£ lj (X) . (2.21) 

fci,..., kn_ x =\ 

Since the coproduct is a morphism, T a (X) also satisfies the defining relations of the Yangian 

Rab{K - Xb) T a {X a ) T b (X b ) = %{X b ) T a {X a ) R ab {X a - X b ) . (2.22) 

Now, we can introduce the main object for the study of spin chains, i.e. the transfer matrix 

M 

t(X)=tr a (T a (X)) = J2Tu(X) . (2.23) 

i=i 

Equation ()2.22j) immediately implies 

[t(X), t(»)]=0 (2.24) 

which will guarantee the integrability of the models (see section 12. 6. 

Let us remark that, at that point, the monodromy and transfer matrices are algebraic objects 
(in (y(gl(Af)))® e ), and, as such, play the role of generating functions for the construction of mon- 
odromy and transfer matrices as they are usually introduced in spin chain models. The latter will 
be constructed from the former using representations of the Yangian, as it will be done below. 

2.4 Symmetry 

The algebraic structure defined above is sufficient to determine the symmetry of the transfer matrix. 
Indeed, we have: 

Proposition 2.2 The gl(N) algebra is a symmetry oft(X). Its generators are expressed in terms of 
the local gl{H) generators as 

T« = L[f ® l^" 1 + 1 ® Lf) ® l^- 2 + • • • + l^" 1 ® L[f . (2.25) 
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Proof: Taking the trace in space a of the exchange relations (|2.22|) . we obtain 



(A« - A 6 ) [ t{X a ) ,T{\ b )) = h{ T(A a ) , T(X b ) ] . (2.26) 
Then, the free term reads 

[t(X a ), T«] = 0, (2.27) 

which proves that = ^fj=i ® -^ij * s a symmetry of the transfer matrix. These generators 
generates the gl(M) Lie algebra. ■ 
Thus, anticipating the spin chain interpretation, we can deduce that all the integrable models con- 
structed in the usual way from t(X) (such as the ones presented in section 12. 7|) possess a gl(M) 
symmetry. In other words, the gl(M) symmetry is valid whatever the Yangian representations are. 
Depending on the model considered (i.e. the choice of representations), we will get the expression of 
the symmetry generators by evaluating the relation ()2.25|) in the representations under consideration. 



2.5 Representations 

As already mentioned, spin chain models will be obtained through the evaluation of the algebraic 
monodromy and transfer matrices in Yangian representations. We thus present here some basic 
results on the classification of finite-dimensional irreducible representations of y(gl(M)). 

2.5.1 Evaluation representations 

Keeping in mind the forthcoming spin chains interpretation, we choose for each local y(gl(M)) 
algebra an irreducible finite-dimensional evaluation representation. 

We start with a finite-dimensional irreducible representation of gl(N)., M(a.), with highest weight 
a = (ax, • • • , ow) an d associated to the highest weight vector v. This highest weight vector obeys 

e kj v = , 1 < k < j < M (2.28) 
e kk v = a k v , l<k<N, (2.29) 

where ai,...,a/v are integers with a k+ i < a k . Indeed the constraints on the parameters a k are 
criteria so that the representation be finite-dimensional and irreducible. Similar criteria will be given 
in Theorem 12.41 for the Yangian. 

The evaluation representation M\(cx) of y(gl(M)) is built from M(ct) and follows from the 
homomorphism ()2.9|) . according to 

L jk (X) v = , l<k<j<H (2.30) 
L kk (\) v = (l - ^j^j v , l</c<Af. (2.31) 

It is important for the following to remark that the previous relations imply that the entries of the 
matrix XC(X) are analytical. 
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The representation M\((l, 0, . . . , 0)), associated to the gl{M) fundamental representation, of C(X) 
provides the R matrix (j2.1|) . 

Let us remark that M\ +a (cx) (a G C) defines also a representation of the Yangian, which is isomorphic 
to M\(et), according to the shift automorphism (j2.17J) . 

2.5.2 Representations of the monodromy matrix 

The evaluation representations of C(X) allow us to build a representation of the monodromy matrix. 
Indeed, evaluating each of the local C a ^ n {\) in a representation M A+aji (a T1 ) for 1 < n < £, the tensor 
product built on 

Mx+^ct 1 )® ■ ■ ■ ® M x+ai (ot £ ) (2.32) 

provides, via ()2.21|) . a finite-dimensional representation for T(A). 

Denoting by v n the highest weight vector associated to ot n = (a", . . . , ajy-), the vector 

v + = v 1 g) . . . ® v e (2.33) 

is the highest weight vector of the representation ()2.32j) i.e. 

T jk (X) v + = , l<k<j<H (2.34) 

/ Pi 71 \ 

T kk (\)v+ = T[ 1-— ^)v + , \<k<M. (2.35) 

n=l V X + a nJ 

We will be interested only in the irreducible finite-dimensional representations of the monodromy 
matrix. When the representation is reducible, the Bethe Ansatz does not give all the eigenvalues of 
the transfer matrix. 

There exists a necessary and sufficient criteria for a tensor product of Yangian representations to 
be irreducible. It uses the following definition 

Definition 2.3 Let X and Y two disjoint finite subsets ofL. X and Y are crossing if there exists 
Xi,x 2 G X and y%,y2 G Y such that 

x 1 < y 1 < x 2 < y 2 or y 1 <xi<y 2 <x 2 . (2.36) 

Otherwise X and Y are non-crossing. 

We associate to each highest weight a the following subset of Z 

X a = {a 1 ,a 2 -l,...,a M -Af+l} . (2.37) 
The theorem giving the criteria to obtain irreducible representations states: 
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Theorem 2.4 [55] 

1. The tensor product M\ +ai (ex 1 ) <g> • • • ® M^ +a£ (a. 1 ) is irreducible if and only if the tensor product 

M^a 1 - %1) ® • • • ® M A (c/ - ^1) 
zs irreducible, with 1 = (1, . . . , 1). 

The tensor product M^a: 1 ) ® • • • <g) M A (ai £ ) is irreducible if and only if all the tensor products 
M\(cx p ) eg) M\(cx q ) with p < q are irreducible. 

3. The tensor product M\(ot) <g> M\(f3) is irreducible if and only if the sets X a \Xp and Xp\X a 
are non-crossing. 

Note that if ex. 1 = a 2 = . . . = ex 1 = ex the tensor product M\(ct)® is irreducible. This special case 
is the one generally used for spin chains models (see examples below). 

2.6 Analytical Bethe Ansatz 

We now use the above mathematical framework to study general closed spin chains. We will be able 
to construct and study a spin chain with arbitrary (and not necessarily identical) representations of 
gl(M) on each site of the chains. Put in other words, the algebraic set-up given above allows us to 
treat simultaneously all the possible spin chain models built in this way. In particular, we will obtain 
the Bethe equations for all these models. 

2.6.1 Hamiltonian of the model 

From now on, we use as local and monodromy matrices the following elements: 

i 

£ a ,n(A) = (A + a n )£ a ,„(A) and f(A) = \[{X + a n ) T(A) . (2.38) 

n=l 

The normalisation of the monodromy matrix ()2.38|) ensures its analyticity. Such a condition is 
crucial for the analytical Bethe Ansatz method. The transfer matrix will be accordingly normalised: 
t(X) = tr a f a {\). 

The properly normalised transfer matrix is a monic polynomial in A of degree t. t(A) = X e + 

e-i 

H n \ n . The £ generalised Hamiltonians H n are in involution (see equation ()2.24|) ) and algebraically 

n=0 

independent (proved by looking at the number of involved sites in each H n ). The Hamiltonian of 
the spin chain model under consideration will be constructed as a polynomial in the generalised 
Hamiltonians H n and will be then integrable. 

Usually, in the spin chain context, we deal with Hamiltonian describing a local interaction, i.e. an 
interaction between nearest neighbour. Unfortunately, at this stage, there is no explicit formula to 
compute this type of Hamiltonian from the transfer matrix. Note however that when all the quantum 
spaces correspond to the same representation, an approach using the fusion of auxiliary spaces can 
be applied [10]. Explicit forms of Hamiltonians will be also given for various models in section l2~7l 
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2.6.2 Highest weight vector / Pseudo-vacuum 

We now compute the eigenvalues of the transfer matrix t(X). As a by-product, they will provide the 
Hamiltonian eigenvalues. The procedure consists in three steps. 

The first step consists in finding a particular eigenvector (so-called pseudo- vacuum) of the transfer 
matrix and in computing the corresponding eigenvalue. We get 

f jk (\) v + = , 1 < k < j < N (2.39) 

i 

f kk {X)v + = Y[(X + a n -ha n k ) v + , 1 < k < AT , (2.40) 

71=1 

where v + is given in (j2.33|) . 

In the following, we use the following notation, for 1 < k < M 

i 

P k (X) = H(X + a n -ha n k ) . (2.41) 

n=l 

These polynomials, called Drinfel'd polynomials, are usually introduced to classify the representa- 
tions of Yangians. 

The highest weight vector ()2.33|) is obviously an eigenvector of the transfer matrix. Indeed, one 
gets 

Af 

t(X) v + = J2 T kk (X) v + = A°(A) v + (2.42) 

k=l 

with 

N 

A°(A) = J2 P k( X ) ■ ( 2 - 43 ) 

k=l 

Note that A°(A) is analytical. In the context of the spin chains, the highest weight vector v + is called 
the pseudo-vacuum. The second step consists in the Ansatz itself which provides all the eigenvalues 
of t(u) from A (A). 

2.6.3 Dressing functions 

We make the following assumption for the structure of all the eigenvalues of t(u) 

A(A) = J2 p kW D k (X) , (2.44) 

k=l 

where D k (X), the so-called dressing functions, have to be determined. At that point, the irreducibility 
of the representation is a necessary criteria for the completeness of the spectrum obtained by dressing. 
From the asymptotic behaviour (A — > +oo) of t(A), we deduce that, for 1 < k < M 

D k {X) ► 1 . (2.45) 
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We suppose that the dressing functions are rational functions of the form 

= n t A V-"" .o-„ n t A t»r ». ■ < 2 - 46 > 

n=1 A A n 2 n=1 A A„ 2 

where M (0) = AfW = 0. 
Remarks 

1. The relation between D^X) and Dfc +1 (A) poles is the basic ingredient of the analytical Bethe 
Ansatz. This pole structure is the simplest one which ensures the analyticity of the eigenvalues. 

2. We introduced shifts in the denominators for later convenience. 

3. The Lie algebra gl(N) being an invariance of the transfer matrix (see proposition 12.2)1 . the 
transfer matrix eigenvectors are indeed eigenvectors of the gl{M) Cartan generators. The 
numbers (1 < k < A/" — 1) are deduced from the action of these Cartan generators on the 
eigenvector of eigenvalue A(A). 

We now tackle the third step, which consists in finding constraints to determine it„ and Vn in terms 
ofA^. 

2.6.4 Fusion procedure 

We shall use the fusion introduced previously in [56,57] to obtain constraints on the dressing functions. 

Let Aj\f be the antisymmetriser defined by the relation (|2.12)) which acts on auxiliary spaces 
di, . . . , aj\f. Then, from the following relation 

f m {\-W+h)---f ai {\) = qdet?(\)A Ar + T a ^\-W + h)---? ai (\)(l-A M ), (2.47) 

we deduce, by taking the trace in the spaces cti, . . . , a^f, that 

t(X - hhf + h) t(X - W + 2h) ... t{\) = qdet f(A) + t f (X) , (2.48) 

where t f (A) = tr ai ^ aAr T a j^(X — hN + H) ■ ■ ■ T ai (X)(l — A^f) is the so-called fused transfer matrix. 



We can compute the value of the quantum determinant using (j2.11j) and the properties of the highest 
weight. Indeed, 

qdet T(A) v + = s 3< a ) Ti,*{i){\ - W + H) • ■ -2aV (A0 (A) v + (2.49) 

N 

= \\ P k (\-W +hk) v + . (2.50) 

k=l 

The quantum determinant being central, the above relation implies that 

AT 

qdet f (A) = Yl ^fc(A -hN + hk). (2.51) 

k=l 
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Then, acting with any eigenvector v with eigenvalue A(A) on relation (|2.48jl . one obtains 

N 

A(A - hj\f +%)... A(A) = \\ P k (X -W+hk) + A f (A) , (2.52) 

k=l 

where Aj(A) v = tf(A) v. Let us remark that this relation shows that v is also an eigenvector of £f(A), 
in accordance with the commutator 

[t f (A),*O)]=0. (2.53) 

Finally, picking the term proportional to Y\^ =1 Pk{X — hM + h~k) in the relation (|2.52jl . we deduce a 
constraint between the dressing functions, namely 

D x {\-W+K) ...Dj<{\) = 1 . (2.54) 

This constraint allows us to express the parameters Un and Vn in terms of An . We conclude that 
the dressing functions take the following form 

M^- 1 ) , _ x(fc-l) _ h(k+l) M< fc > , .(*) h (fc-2) 

= n n : >; a • ^ 

n=l A A n 2 n=1 A A n g 

2.6.5 Universal Bethe equations 

We have chosen the normalisation of the matrix T(A) in such a way that its entries are analytical. 
Then, the eigenvalues of t(A) are also analytical, since t(X) can be diagonalised by a constant matrix 
(see equation ()2.24|) ). 

Theorem 2.5 The Bethe equations read, for 1 < k < JV — 1 and 1 < n < M^) 

M^" 1 ) M« M( fe+1 ) p I _|_ 

n e - w fc) - n w fc) - ^ ] ) n ^ - a ^ +i) ) = 2 ^ ( 2 - 56 ) 

m=l m=i m=l ( "I - n 



where 



e x (X) = • (2-57) 

The left hand side of \2. 56]) depends only on the choice of the algebra (the indices of the function 
e x (X) describe the entries of the Cartan matrix of gl(N)), while the right hand side depends on the 
choice of the representation. 

Proof: By imposing that the A(A) residue vanishes at A = An + ^r, we find (12.56(1 . ■ 
The RHS of (|2.56() can be written in terms of the functions e x (X), using the expression of the 
highest weights. These Bethe equations have been computed in [12], however the method and the 
starting hypotheses are different. The identity between the results appears as a ground for this 
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Ansatz. 



It should be clear that the Bethe equations ([2.560 . and the dressing of the eigenvalues, ([2.44(1 and 
([2.550 . are valid whatever the expression of the Drinfel'd polynomial is, and as such are universal. 
The dressing functions (and thus the expression of the eigenvalues) appear formally independent from 
the choice of the representations. However, the Bethe equations depending on the representations, 
their resolution will lead to different eigenvalues. 

The choice of a closed spin chain model amounts to the choice of the gl(M) representation M(a k ) 
for spins at sites k, 1 < k < £. This will fix the evaluation representations M\(ct k ), hence the 
polynomials -Pfc(A). Then, the eigenvalues and the Bethe equations follow. We now illustrate this 
procedure by employing spin chain models. 

Remark: reducible representations 

When the representation is reducible, the above calculations are still valid, but they do not lead to a 
complete set of eigenvalues for the transfer matrix. In fact, one gets in that case all the eigenvalues 
associated to the irreducible subrepresentation with highest weight v + . A simple indication for that 
is the eq. ([2.50)1 which now implies ([2.51(1 only on this irreducible subrepresentation. 

2.7 Examples 

Choosing appropriate representations, we shall recover known results associated with the fundamental 
representation, generalise the relations about the alternating spin chains and provide new integrable 
models (such as general impurity spin chains). For simplicity, we will most of the time set the 
inhomogeneous parameters a n to zero. However, our formalism easily deals with these inhomogeneous 
parameters, as we shall see in the next example. 

2.7.1 Closed spin chain in the fundamental representation 

The usual closed spin chain corresponds to spins in the fundamental representation. The Hamiltonian 
is given by the well-known formula 



In this case, we have a n — (1, 0, . . . , 0), for 1 < n < i. Then, the Drinfel'd polynomials read 

( I 



H 



d 



(In t(X)) 



A=0 



(2.58) 



Y[(\ + aj -h) , k = i 



= { 



3=1 
I 



(2.59) 



so that 




(2.60) 
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Plugging these expressions in the Bethe equations (|2.56j) . we recover the usual Bethe equations for 
closed spin chains 3 



n e-i (a? ) - as-^) n ^ w* } - ^ n -w-e^ M e ' 

m=l m=i m=l I -, r / 1 

Since the value of the local operator £y(A) at A = is the permutation operator between spaces 
i and j (see eq. (|2.2J) ). we can construct a local Hamiltonian by the relation (|2.58|) (when a n = 0) 

i 

H<xJ2 V n~l,n with Vqi = V ei . (2.61) 

n=l 

In the case of gl(2), one recovers the celebrated XXX Hamiltonian 

1 - 

H = 2 E ( + CT n-i^ + <-i< + 1 ) , (2-62) 

n=l 

where cr* , t — x,y,z are the Pauli matrices at site n, and ctq = o\. 

If we take a p ^ for a particular site p (and a n = for n ^ p), we obtain a Hamiltonian with 
one impurity 

l 

> ^ _ f}/ 1/2 , . . 'i ) \ /^n t"> \ 

if OC 2^ K-l,n - - _ ^ + a 2 _ft2 \ a p^P-hP+l - h r p+ i,p ) + q2 _ ^P-hP+n^P-hP ~ rp+l,p) ■ 

n—l P p p 

2.7.2 Closed spin chain for non-fundamental representations 

One can generalise the above example to the case where all the spins belong to the same (not 
necessarily fundamental) representation, given by 

a 1 = a 2 = . . . = a £ = (cti, a%, . . . , a^i) • (2.63) 

In particular, we recover the result given in [10,11,23,58,59] about the XXX higher spin chains. 
We will use the variables 

a k ± a k+1 = $ 

which are integers, since we consider gl{M) irreducible finite-dimensional representations. This leads 
to the following Drinfel'd polynomials 



P k (\) = (A - ha k ) e so that v 



P l \^ _L hk 
-ffc+l I A n + — 



e«- f Al fe ) + ft* ~ K 



;2.64) 



3 We remind that h = —i when dealing with spin chain models. 
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For instance, if we particularise to the gl(2) spin chain in the spin s representation, we get as Bethe 
equations 



m r , 1 —2s\ 

Y\ e 2 (K - A m ) = e 2s ( X n + h — - — J 

m — i L V / 



M 

(2.65) 



The construction of a local Hamiltonian cannot be repeated from section 12 . 7. II because there is no 
particular parameter where the local operator C(X) is the permutation. However, a local Hamiltonian 
can be constructed by using the fusion method introduced in [10,11,23,58] or by evaluating the 
universal .R-matrix, see e.g. [59]. It takes the form 

l 2s / j \ 2s _ 

H oc Qis ( Ci< + + ) , where Q 2s (x) = J] E T II ■ ( 2 - 66 ) 

n=l j=l \fc=l / «=o J ' 

In the above formula, X[ = ^(1(1 + 1) — 2s(s + 1)), s^, t = x,y, z are the gl(2) generators in the spin 
s representation acting in the quantum space n, and satisfying Sq = s\. The energy spectrum is then 
given by 

j=i j 

where Xj are solutions of the Bethe equations (|2.65|) . 
2.7.3 Alternating spin chains 

In alternating spin chains, the spins along the chain belong alternatively to two different given 
representations. We take the particular example of the alternating spin chain with the number of 
sites i = 2i even. The spins of even sites are represented in the fundamental representation whereas 
the spins of the odd sites are in another representation. We take the following particular example 
where the highest weights are given by 

a n _ \ (1,0, • • • ,0), \<n<t and n even 
° [(2,0,...,0), \<n<i and n odd 

Then, the left hand side of the Bethe equations read, for 1 < n < , 1 < k < Af — 1 and J\f > 2 



ei ^je 2 ^-fjj , k = l {2m) 
P k+ i[X { n k) + ^ [i , Kk<Af 

We recover the Bethe equation given in [8] (see also [60]) for gl{2). In the case gl(2), we can compute 
a Hamiltonian by the usual formula ()2.58|) which contains both nearest and next-to-nearest neighbour 
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interactions with periodic boundary conditions. The explicit form of this Hamiltonian is given by [8] 



H oc (2cr 2j - • s 2j+ i + l)(2cr 2j+2 • s 2j+ i + 3) 

+ (2 1) [(s 2i _a • s 2j+1 + l)(2<r 2j ■ a 2j+1 + 1) + 2] } (2.70) 

where a = (a x ,a y ,a z ) are the Pauli matrices and s = (s x ,s v ,s z ) are the generators of sl(2) in the 
spin 1 representation. 

We can also recover the results of [9] where another type of alternating spin chains has been 
studied for su(3). 

2.7.4 Impurity 

We consider now a spin chain with one site (the impurity) in a representation different from the 
others. Let us take as example a spin chain where all sites are represented in the fundamental 
representation except for the p th which is associated to the representation of highest weight a p . In 
this case, the left hand side of ()2.56|) becomes 



(*0 j_ hk, 
2 



P ( \^ _L Hk 
r k+1 I An H 



(2-71) 



\W I H k fc P 

I. An + - 2 " 



fe ^ 1 



The Hamiltonian can be written as 

l 

HCX V n-l,n + ( Vp-l,p+l Cp+lA Q ) ~ h ) £p-l, P (°) ( 2 - 72 ) 



n=l 



where Poi = ^i, ^„ iP (A) = J^ij^ij <8> (A — /i^ji). In the last formula, belongs to the space n 
(fundamental representation), while (A — hSji) is in the particular space p where the generators of 
gl{M)i Sij, are in the representation with the highest weight ct p . 

2.7.5 Generalisation to tensor products of representation on each site 

Up to now, we have assumed that on each site of the spin chain, only one evaluation representation 
occurs. This assumption is natural from the spin chain point of view, since one can interpret the 
underlying gl{M) representation as carrying the spin. However, the algebraic framework we have 
presented can deal with more general representations of the Yangians, provided they are irreducible 
and finite-dimensional. Using the theorem 12 A\ the irreducible representations will be constructed 
from tensor products of evaluation representations. Let us stress that, generically, this tensor product 
of evaluation representations is irreducible , although for the transfer matrix symmetry algebra gl(N) 
these representations are fully reducible. 
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From the physical point of view, the model will describe a spin chain possessing on each site a 
quantum space which is a tensor product of evaluation representations of the Yangian. However, this 
model (in particular the transfer matrix) can be reinterpreted as a usual spin chain model but with 
a higher number of sites, each of them associated to only one evaluation representation. 

Finally, let us remark that this construction is in essence opposite to the fusion procedure. Indeed, 
for the fusion, one takes particular points (described in theorem I2.4JI where the tensor product of 
evaluation representations is reducible . 

3 Open spin chains with preserving boundary conditions 

In this section, we compute, along the lines described in the previous section, the Bethe equations 
for the open spin chains with soliton preserving boundary conditions [61-64]. For such a purpose, 
we first need to introduce some new algebraic objects such as the reflection algebra or the K matrix. 

3.1 Reflection K matrix 

In the case of soliton preserving boundary conditions, we need to introduce numerical matrices, called 
K matrices, which are solutions of the reflection (boundary Yang-Baxter) equations [27]: 



R ab (X a - \ b ) K a (\ a ) R ba (\ a + Aft) K b (\ b ) = K b (X b ) R ab {\ a + A 6 ) K a (X a ) R ba (X a - A 6 ) . (3.1) 



The K matrix is interpreted as the reflection of a soliton on the boundary, coming back as a soliton. 
The solutions of the equation (|3.1|) have been classified in [65]: 



Proposition 3.1 Any invertible solution of the soliton preserving reflection equation h8. 1\) takes the 
form K(X) = U ( E + t-^aM U" 1 where either 



(i) E is diagonal and E 2 = 1^- (diagonalisable solutions) 

(ii) E is strictly triangular andW? = (non- diagonalisable solutions) 

The matrix U is an element of the group GL(N) and £ a free parameter. The classification is done 
up to multiplication by a function of the spectral parameter. 

Note that all the K matrices (but zero) obey a relation K(X)K(—X) = f(X) Ij\f for some non-zero 
even function f . 

A suitable relabelling of the indices allows us to choose the matrix E in (i) of proposition 13.11 as 



E = diag(l, . . . , 1 




-1,...,- 




(3.2) 



M 



N-M 



with < A4 < N ' . In the following we will only deal with diagonal solutions of the form 



K{\) = diag( \ + £,. ..,X + { , r A + £,■■■, -A + £ ) . 



(3.3) 



M M-M 

We normalise the K matrix so that its entries be analytical. 



16 



Case of non-diagonal reflection matrices: The general treatment of non-diagonal reflection 
matrices is yet an open problem. In the case where each spin is represented in the fundamental 
representation, the problem has been solved in [56,57] for K + = 1 and in [66] for simultaneously 
diagonalisable reflection matrices K + and K~ . In the case of the XX Z model, the procedure given 
in [56,57] to treat the non-diagonal reflection matrices does not work. However, interesting devel- 
opments have been done in [67,68] to attempt a general treatment of non-diagonal reflection matrices. 



3.2 Reflection algebra 

The reflection algebras are constructed as subalgebras of a Yangian, which is here y(gl(N)). Starting 
from the generators T(A) of y(gl(Af)) introduced in ()2.8|) . we define 

B(\)=T(\)K(\)T(-\y 1 . (3.4) 

£>(A) generates an algebra, denoted %$(J\f, M), whose exchange relations are given by 

R ab (\a - Xb) Ba{K) R b a(X a + A 6 ) B b {X b ) = B b {\ b ) R ab {\ a + A 6 ) B a {\ a ) R ba {\ a - A 6 ) . (3.5) 

Writing £>(A) as 

M +°° »(n) 

B(A) = X)^®S«(A) = X;^r 

i,j=l n=0 

one can show that B^ generates a gl{M) © gl{M — M.) subalgebra in 23 (TV, M). 

Another reflection matrix K + (X), solution of an equation dual to (j3.1j) . is usually introduced to study 

open spin chains [28]. For simplicity, we will take here K + (X) = Iby. 

The (algebraic) monodromy matrix used to construct open spin chain is obtained from the local 
operators C a j{X) of the Yangian (|2.20|) . It takes the following form 

B a (X) = C al (X) . . . O(A) K a (X) C'Ji-X) . . . C-^-X) . (3.6) 

The transfer matrix becomes 

AT 

b(X)=tr a (B a (X)) = J2Bu(X) (3.7) 

i=i 

and, as in the Yangian case, the commutation relations defining the algebra allow us to show 

[6(A),6(//)] = 0. (3.8) 

This relation (|3.8|) guarantees the integrability of the model, usually described by the following 
Hamiltonian 

ff=-~6(A) • (3.9) 
2 a A A=o 

Anticipating again the physical spin chain interpretation, one can compute the symmetry of these 
models 
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Proposition 3.2 The transfer matrix 6(A) describing open spin chain models admits an gl(Ai) © 
gl(J\f — M) symmetry. 

Proof: Following the steps given for the closed spin chains (see proof of proposition 12. 2|) . one shows 
that [£>^ , 6(A)] = 0. Since £>^ generates a gl(M) © gl(N — M) algebra, this ends the proof. ■ 



3.3 Representations of 05 (J\f, M) 
3.3.1 Representation of T _1 (A) 

In order to study the representations of 05 (J\f, M), we start from the representations of the Yangian 
introduced in the section |2~51 Let M\(ct) be an evaluation representation of £(A) with the highest 
weight vector v. We can show that v is also a highest vector of £ -1 (A) with 

L' jk (X)v = , l<k<j<N (3.10) 
L kk\X) v = A tt — r tt — — - — r v , 1 < k < N , (3.11) 

where L'- k {\) are the matrix elements of £ -1 (A). The values appearing in (J3.ll)) are computed from 
the relation (j2~TK|l . 

The relations (J3.10)) and (J3.ll)) imply that v + as given in (J2.33)) is the highest weight vector of T _1 (A) 
with 

T; k (X)v + = , l<k<j<N (3.12) 

T ^ (A) " =A P l( A)---P,(A + ^-^) V ' ( 3 - 13 ) 

where T'- k {\) are the matrix elements of the matrix T -1 (A) and -Pfc(A) are defined in 1)2.41)1 . 
Let us remark that 



Px(-A) P 2 (-A + H) ■ ■ ■ Pm{-\ + W-K) rr _ Xt 

:-av 



T(A) = ^^>*^-^ 7 ^y-^^-"> T -\-X) (3.14) 



can be understood as the Yangian generators represented in the following tensor product of evaluation 
representations (as defined in section 12". 5. 2)1 

i 

(g) (MaOS 1 '") © • • • © M x (f3^ n )) (3-15) 

n=l 

where (3 k,n = (1 — a™, . . . , k — 1 — a k _ v 0,k — a k+v . . . , TV — 1 — ajy). This shows that the matrix 
(_A)(^-iKf(A) is analytical. 
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3.3.2 Representation of the monodromy matrix S(A) 

We can now describe the representations of the monodromy matrix B(X) defined by (jH.ojl . It is 
known [69] that any finite-dimensional representation of 53 (J\f, M) is a highest weight representation. 
They can be constructed in the following way 

Theorem 3.3 Let us consider the Yangian highest weight representation M\(a}) ® . . . (g) M\(a l ) 
with highest weight vector v + = f 1 ®. . .<S)v . Then the realisation \3. 6]) generates a 33(JV, A4) highest 
weight representation, whose highest weight vector is also v + with 

B jk (\)v + = 0, l<k<j<N (3.16) 
B kk {\)v + = %( A ) ^( A ) + 2A -kh + h ^ (A) ) "' + ' 1 < k <^^ ( 3 - 17 ) 



where, for 1 < j < Af, 

(2\-jh+h)(2\-jH) 

and 



-2A 

a jW = , n , , fe ; fcUo , ,^ (3-18) 



/•ifc(A) 



Proof: A direct calculation (similar to the one done in [69]) leads to, for 1 < k < A4, 
2\-kh+h„ h '~ 



(3.19) 



B kk (\)v + + — %( A K = (A + T kk {\)T' kk {-\)v + (3.20) 



2A y ' 2A 

and for M + 1 < k < M 

o\ — kh -\- ft ft 

2A + B kk (X)v + + £-Y. B ^ X > + = (" A + * + ^ T ^( A ) T ^(- A )^ + • ( 3 -21) 

Then, inverting these formulae and using the expressions ()2.35|) and ()3.13|) . one gets the expression 
3.4 Analytical Bethe Ansatz 

The analytical Bethe Ansatz method is based upon the analyticity of the represented generators of 
the algebra. It is ensured by a suitable normalisation given in the following proposition 
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Proposition 3.4 Let 

g(X) = (-1) £ p^-X^.-.p^-X+Mh-h) B(X). (3.22) 

Then, B(X) is analytical (in X). 
Proof: B(X) can be rewritten as 

B(X)=f(X) x K(X) x ((-A) ( -^- 1)£ T(A)) . 

The three terms of this product are analytical. ■ 
From now on, we will use £>(A) instead of B(X) to ensure, as in the closed spin chain case, the 
analyticity of the eigenvalues of the transfer matrix 



b(X)=tr a B a (X). (3.23) 



3.4.1 Pseudo- vacuum 



As in the case of the closed spin chain, the first step of the analytical Bethe Ansatz consists in finding 
a particular eigenvalue of the transfer matrix. This eigenvalue is computed thanks to the highest 
weight vector v + . Indeed, one gets 



b{X) v + = J2 B kk {\)v + = A°(A) v + (3.24) 
fc=i 

where 

A°(X) = J29k(\)f3 k (\). (3.25) 

fc=i 

The functions #fc(A) depends only on the boundary matrix. They are given by 

2A(2A-A^) K ( A + e f ° r l ^ k ^ M (3 26) 

(2A /,•// • /;)(2A /,•/;) X | A • i • Mh for M + l<k<N. 

The functions f3 k (X) depend on the choice of the representation: 

(3 k {X) = p 1 (-x + h)--- P fe _i(-A + kh- h)p k (X)p k+1 (-x + kh)--- Pm(-x +Nh-h). (3.27) 

Let us remark that A (A) is analytical and in particular its residue for A = kh/2 vanishes (0 < k < Af). 
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3.4.2 Dressing functions 

The central hypothesis of the analytical Bethe Ansatz is that all the eigenvalues of 6(A) can be 
written 

AT 

A(\) = ^ / g k (\) /3 k (\) D k (X), (3.28) 

k=l 

where the dressing functions D k (\) are rational functions and need to be determined while g k (X) and 
/3fc(A) are given by ()3.2(ij) and (|3.27)l . respectively. The vanishing of the residues of A(A) at A = kh/2 
implies that 

D k (kh/2) = D k+1 (kH/2) for 1 < k < M - 1 . (3.29) 

Starting from the expression (|2.4fi|) for the dressing functions, one can show that the M^'s are even, 
and that (up to a rescaling — > /2) the dressing functions read 



D k (\) 



n 

n=l 


A 


+ A^ 


fi(fc+i) 

2 


A 


x(fc-l) R(fc+1) 
An 2 


A 


+ xr i] 


ft (Jfc-1) 
2 


A 


R(Jfe-l) 
An 2 


MO) 
n=l 


A 


+ A n — 


f + ft 


A- 


A«-f + £ 




A + Ai fc) 


hk 
2 


A 


-A„ - — 



(3.30) 



where M (0) = = 0. 

3.4.3 Bethe equations 

The normalisation of the matrix i3(A) has been chosen in such a way that its entries are analytical. 
Then, the eigenvalues of the transfer matrix 6(A) are also analytical (since the diagonalisation matrix 
does not depend on A). 

Theorem 3.5 The Bethe equations read, for 1 < k < M — 1 and 1 < n < 

n M^ur^n^A^ a«) n m^u^) 

m=l m =! 771=1 

Pk(\ i n ) + ! f) f-e-M-v/n^) if k = M 

x \ (3.31) 



Pk+i ^Ai ^ + y\ otherwise 

where 

e x (X, fi) = e x (A - /i) (A + //) , (3.32) 
the functions e x (A) are defined by flO^ and M<°) = AfW = 0. 

Proof: By imposing the vanishing of the A(A) residue at A = An + Ap, we obtain (J3.31|) . ■ 
As in the case of the Yangian, the left hand side of ()3.31|) depends only on the choice of the algebra 
whereas the right hand side depends on the choice of the representation and the K matrix. 
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3.5 Examples 
3.5.1 Generalities 

All the cases presented in section 12.71 can be treated in a similar way for the open spin chain, using 
the usual formula given in [28] for the Hamiltonian. 

As a basic example, one can easily check that the present approach reproduces correctly the 
results obtained for the open gl(J\f)-spm chain with generic boundary [57]. 

As more involved examples, we can generalise directly the spin s chain and the alternating spin 
chain (see sections [2.7.21 and 12. 7. 3j) by adding a boundary with the procedure given above, extending 
the results obtained in [70]. 



3.5.2 Boundaries with operators 

One may wonder whether the boundary matrices K± can be promoted to operators. Indeed this 
amounts to "fuse" the boundary to the last site to get a dynamical boundary. This was considered 
for instance in [28,71-74]. 

We treat here an example suggested by K. Zarembo. We study the gl{2) spin chain with t — 2 
spins 1 in the bulk and for two spins 1/2 the boundaries. For this spin chain, we represent the 
monodromy matrix where the highest weights are given by 

Q »=/ (1 ' 0) • B = 1 '< (3.33) 

I (2,0) , l<n<£ v 1 

The corresponding integrable Hamiltonian is given by (up to an irrelevant overall normalisation) 



1-2 , v 

H = 2(7, ■ S 2 + 2S L -x ■ <r L + ]T [Si ■ S i+1 - - (Si ■ S l+1 ) 2 j 



(3.34) 



i=2 

with the following conventions: 



"-O. ^-fiS). -"ft- -!. ("0 



/0 1 0\ /0 0\ (2 \ 

S + = 1 , S~ = V2 1 , S z = io , (3.36) 
\0 0/ \0 1 0/ \0 -2/ 

and 

A B = 2A + B- + 2A"fi + + A Z B Z . (3.37) 

The S used here is actually equal to 2s of section 12.71 so that cr and S have the same commutation 
relations. 

This Hamiltonian comes from the monodromy matrix ()3.4j) with the following prescriptions: the 
auxiliary space is three-dimensional (spin 1 representation of gl(2))\ the reflection matrices are 
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taken to be identity matrices. 

Hence 6(A) = Tr a T a (A)T a (-A)- 1 with 



T(A) = ^' i) (A)</- ) i(A) • • -R^W^iX) (3.38) 



and 



R%*\\) = AI3 ® I 2 — ^S a ■ a-j j — 1,£ , (3.39) 
<' 1} (A) = {X + h) ^~ 2h) ^h-^S a .S j + ^(S a . Sj r , J = 2,... ,£-1.(3.40) 

These i?-matrices can be derived from -RjJ' 2 given by ()2. 1|> using the usual fusion procedure [75]. 
As a consequence, the transfer matrix 6(A) (and the Hamiltonian) commutes with the transfer matrix 
6(A) built with the same quantum spaces and spin 1/2 auxiliary space. 

Then, using (|3.31|) . the Bethe equations are, for 1 < n < M 

M 

I [ e 2 (A„ - A m ) e 2 (A n + A m ) = e x (A n ) £+2 e 3 (A n /~ 2 (3.41) 

m— 1 

The bulk part of 1)3.34)1 is the mixing matrix for some sort of gluon operators in large-A QCD, 
see [45]. The spin chain boundary term in (|3.34|) corresponds to the quark-gluon operators. 

Perspectives 

A natural development of this work is the generalisation to the trigonometric case, which will be 
presented in a further publication. Soliton non-preserving boundary conditions will also be studied 
in this framework. 
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